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Context. One of the often discussed models for high-frequency quasi-periodic oscillations of X-ray binaries is the oscil- 
lating torus model, which considers oscillation modes of slender accretion tori. 

Aims. Here, we aim at developing this model by considering the observable signature of an optically thick slender ac- 
cretion torus subject to simple periodic deformations. 

Methods. We compute light curves and power spectra of a slender accretion torus subject to simple periodic deforma- 
tions: vertical or radial translation, rotation, expansion, and shear. 

Results. We show that different types of deformations lead to very different Fourier power spectra and therefore could 
be observationally distinguished. 

Conclusions. This work is a first step in a longer term study of the observable characteristics of the oscillating torus 
model. It gives promising perspectives on the possibility of constraining this model by studying the observed power 
spectra of quasi-periodic oscillations. 
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1. Introduction 

Some microquasars exhibit high-frequency (millisecond) 
quasi-periodic oscillations (QPOs), which in this work al- 
ways refer to high-frequency oscillations; low-frequency 
QPOs are not considered here. These high-frequency oscil- 
lations are characteriz ed by a narrow peak in the source 



power spectrum (see van der 



Klis 
~Ks~t 



2004 Remillard & 



|McClintock| |2006| for a review) 

variability is of the same order as the Keplerian orbital 
period at the innermost stable circular orbit of the cen- 
tral black hole, it is highly probable that these phenomena 
are linked with strong-field general relativistic effects. The 
other different low-mass X-ray binaries that contain neu- 
tron stars display QPOs on these prominent timescales as 
well. There is no consensus on the possible QPO origin and 
uniformity across different sources yet. 

A variety of models have been developed so far to ac- 
count for t he high-frequency variability. |Stella fc Vietri| 
( 1999[ ) and |Stella et aL| ( |1999| ) suggest that QPOs could 
arise from the modulation of the X-ray flux by the pe- 
riastron precession and the Keplerian frequency of blobs 



of matter orbiting in an accretion disk around the cen- 
tral compact object. Also QPOs could be due to modu- 
lation of the X-ray flux by oscillations of a thin accretion 
disk surrounding the central compact object (see |Wagoner| 
[19991 pCatc][200l| . |Fragile et ah] pOOlj ) propose that QPOs 
are due to the modulation of the X-ray flux caused by a 
warped accretion disk surrounding the central compact ob- 
ject. Poi nting out the 3:2 ratio of some QPOs in different 



e timescale of this sources, Abramowicz & Kluzniak ( |2001[ ) have proposed a 



resonance model in which these pairs of QPOs are due to 
the beat between the Keplerian and epicyclic frequencies of 
a particle orbiting around the central compact object. 

General relativistic ray tracing of radiating or irradiated 
hot spots orbiting around black holes and neutron stars 
was founded at the turn of the millennium. For instance, 



Karas ( 1999 ) discusses already sophisticated numeric mod- 
eling of observable modulation from clumps distributed 



around certain preferred circular orbits. Later, Schnittman 



& Bertschinger (2004) investigated in great detail predic- 
tions of a model of hot spot radiating isotropically on nearly 
circular equatorial orbits. In their study, the hypothetic res- 
onance between the Keplerian and radial epicyclic frequen- 



1 



Mazur et al.: Oscillating tori for modeling QPOs 



cies gives rise to peaks i n the modeled power spectrum. 
Tagger fe Varniere] ( 2006 ) advocate the fact that QPOs in 



microquasars are due to the triggering of a Rossby wave 
instability in the accretion disk surrounding the central 
compact object. Ray-traced light curves have been recently 
developed for this model by Vincent et al. (2012 1. Finally, 
models involving oscillations of accretion tori have been de- 
veloped. As our work is a follow-up of these past studies, 
we present them in more details. 

The first study of a QPO model invo lving thick ac- 
cretion structure (tori) was developed by Rezzolla et al 
( 2003 ), who showed that p-mode oscillation of a numerically 



computed accretion torus can generate QPOs. Ray-traced 
lig ht curves and power spectra o f this model were derived 



by Schnittman & Rezzolla 
accretion tori (for which 



fl2006j . The analysis of analytical 
physical quantities are known 
analytically throughout the struct ure and at al l times ) as a 
model for QPOs was initiated by Bursa et al. ( 2004 ), who 
performed simulations of ray-traced light curves and power 
spectra of an optically thin oscillating slender torus. This 
model took only into account simple vertical and radial si- 
nusoidal motion of a circular cross-section torus. To allow a 
more general treatment, a series of theoretical works were 
dedicated to developing a proper model of general oscilla- 
tion modes of a slender or non-slender perfect-fluid hydro- 



dynamical accretion torus (Abramowicz et al. 2006 Blaes 



et al.||2006| |Straub fc Sramkova|2009| ) . Considering hydro 



dynamical perfect-fluid tori is interesting in terms of sim- 
plicity as all computations can be derived analytically and 
all physical quantities are known analytically at all times. 

The aim of this article is to strengthen the perfect- 
fluid oscillating slender torus model for QPOs by progress- 
ing towards determining its observable signature. Here, we 
present a first step towards this goal, which consists of 
determining the observable signature of simple deforma- 
tions of a slender torus in the Schwarzschild metric. We 
do not use the ma thematically well-d efined oscillation for- 
malism derived by Blaes et al. (2006). We rather consider 
the impact on the observables of elementary deformations 
of the torus cross-section: translation, rotation, expansion, 
and shear. We simulate light curves and power spectra of 
such tori subject to these simple deformations, taking into 
account all relativistic effects that will affect radiation by 
using a ray-tracing code. The focus of this article is thus to 
determine the impact of the geometrical change of a slender 
torus on the observed flux variation. We do not consider the 
physics of the matter forming the torus. 

This model is extremely simple and does not claim to 
give a realistic view of the nature of QPOs. However, it does 
allow determination in the simplest possible framework of 
whether different kinds of motions of a toroidal accretion 
structure lead to potentially observable differences in the 
light curves. This work thus makes it possible to go one step 
further in the development of the oscillating torus model for 
QPOs. 

Section [2] describes the equilibrium slender torus, while 
Section [3] describes its deformations in terms of translation, 
rotation, expansion, and shear. Section [4] shows the ray- 
traced light curves and power spectra of the deformed torus, 
and Section [5] gives conclusions and prospects. 



2. Equilibrium of a slender accretion torus 

This section derives the equations that describe the accre- 
tion torus at equilibrium. The spacetime is described by 
the Schwarzschild metric in the Schwarzschild coordinates 
(t, r, 6, <p), with geometrical units c = 1 = G and signature 
( — ,+ ,+,+). The line element has then the following form: 



ds 2 = gttdt 2 + g rr dr 2 + ; 



2M\ , , 
I"— J d* 2 + 

; (d9 2 + sin 2 6»d^ 2 



d6> 2 
1- 



i) 



where M is the black hole mass. In the remaining article, 
we use units where this mass is 1. 

In this spacetime, we consider an axisymmetric, non- 
self-gravitating, perfect fluid, constant specific angular mo- 
mentum, which circularly orbits the accretion torus. This 
torus is assumed to be slender, meaning that its cross sec- 
tion is small as compared to its central radius. 

2.1. Fluid motion 

As the fluid follows circular geodesies, its 4- velocity can be 
written as 

u" = A(r]i" + n^), (2) 

where f2 is the fluid's angular velocity, and -rf and £ M are the 
Killing vectors associated with stationarity and axisymme- 
try, respectively. The constant A is given by imposing the 
normalization of the 4- velocity, u^u^ = — 1. 

As spacetime is stationary and axisymmetric, the spe- 
cific energy £ and specific angular momentum C are defined 
as 

£ = -rfup = -u t (3) 

which are geodesic constants of motion. 

The rescaled specific angular momentum £ is defined 
according to 

e=~ = -^-. (4) 

£ u t 

We assume this rescaled specific angular momentum to 
be constant throughout the torus: 



const. 



(5) 



The 4-acceleration along a given circular geodesic fol- 
lowed by the fluid is 



1LA 



(0) 



where U is the effective potential (see, e.g., Abramowicz 
eTaI][2006l ): 



U = g tt +i 2 g 



(7) 



The radial and vertical epicyclic frequencies for circular 
motion are related to the second derivatives of this poten- 
tial: 

1 £ 2 



w 2 e 



' d 2 u 

2 A 2 g rr V dr 2 

1 £ 2 ( d 2 u 

2A 2 g ee \d6 2 



(8) 
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In the Schwarzchild metric, these epicyclic frequencies 
at the torus center are 



J r0 



O 2 



6 
''0 



-no 



= CI 



K 



(9) 



where a subscript denotes here and in the rest of this arti- 
cle a quantity evaluated at the torus center. The Keplerian 

3/2 

angular velocity is well known: VLk = r . 

In the rest of this article, the torus central radius will 
be fixed to the value ro, satisfying the following condition 



3.1. Deforming the torus 

In this section, we consider various simple time-periodic 
deformations of the torus cross-section. Our aim is to de- 
termine the new torus surface function corresponding to 
these deformed states. 

All deformations boil down to performing transforma- 
tions on the (x,y) coordinates of the form: 



x(t) 

m 



ai(t)x{0) 



■a 2 (t)y(0) + a 3 {t), 
b2(t)y(0) + b 3 (t), 



(18) 



^00 
L0 r o 



(10) 



This choice is linked to our goal of applying the de- 
formed torus model to twin-peak QPOs. It results in 



r = 10.8. 



(11) 



2.2. Surface function 



Following Abramowicz et al. (2006), we define the surface 
of the torus by the locus of the zeros of a surface function 
f(r,6). Introducing the following set of coordinates 



r - r 
ro 



V = (Vgoe)r 



tt/2 - 
ro 



(12) 



Abramowicz et al. ( 2006 1 show that the surface function is 
expressed according to 



where 



/ = 1 



U! r o 



1 



1-2 2,-2 2\ 
K X +UJ 90 y ), 



UJ r o 



ueo 



Woo 



(13) 



(14) 



the parameter /3 is related to the torus thickness, and the 
torus being slender: 



/?< 1. 



(15) 



Abramowicz et al. (20061 also show that the x and y 
coordinates are of order fj. We thus define a new set of 
order-unity coordinates: 



x 



and: 



f = l- (u, 2 x 2 +u, 2 y 2 ) 



(16) 



(17) 



The equilibrium slender torus has therefore an elliptical 
cross section. 



where the and 6; functions are simple trigonometric func- 
tions of time. Table [l] gives these functions for all defor- 
mations considered in this article as a function of a free 
parameter A and of the deformation's pulsation u). This 
pulsation will have the same value for all the simulations 
performed in this article. We choose to = CIk- We note that 
this value is a choice as the deformations we consider are 
not physically justified, but imposed for the simplicity of 
the model. We make here the simplest and most natural 
choice of considering the Keplerian pulsation. 

At this stage, the deformed torus we are considering has 
thus only two degrees of freedom: 

- the torus thickness j3 parameter, with (5 <C 1, 

- the deformation parameter A, with typically A ~ 1. 

The actual values of (3 and A will be chosen to satisfy 
the slender torus approximation for all times. 

3.2. Perturbed fluid 4-velocity 

The fluid 4-velocity in the deformed torus is given by 



,n — 



(19) 



where Uq is the equilibrium 4-velocity already defined in 
Eq. (pi), which is assumed to be everywhere equal to its 
central value, in the slender torus limit. 

The value of the perturbation 5u^ induced by the defor- 
mation is easily computed by using the coordinate trans- 
formations given in the previous Section. By writing a first- 
order expansion of x(t + 5t) and y(t + St) for a small time 
increment St, it is straightforward to obtain the expressions 
of dx/dt and dy/dt, which are themselves linearly related 
to dr/dt and dO/dt. 

An example of this is the case of radial translation of 
the torus cross-section: 



u r = Su r 

Ug = Sug = 0, 



(20) 



3. Deformations of a slender accretion torus 

This section is devoted to determining the surface function 
of the torus subject to simple deformations: translation, 
rotation, expansion, and shear. The 4-velocity of the per- 
turbed fluid is also given. 



with 



ou r = = 
dt 



XiOTQ 
,V 9rr)r 



cos(o;t)u t 



(21) 



from Table [U 

As the torus is assumed to have a constant rescaled an- 
gular momentum £q> the following relation holds between 
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Deformation ai <22 0,3 b\ 62 63 

Radial translation 1 OA sm{ut) 1 

Vertical translation 1 1 A sm(ujt) 

Rotation cos(tot) sin(urf) — sin(ujt) cos(ujt) 

Expansion 1 + Asin (u;t) 1 + Asin(aii) 

Radial simple shear 1 Asin(a;i) 1 

Vertical simple shear 1 Asin(ut) 1 

Pure shear 1 + Asin(wt) [1 + Asin(^)]" 1 



Table 1. Functions and hi of Eq. 18 for all periodic deformations, A is a free parameter and u> is the deformation's 
pulsation. 



the non-equilibrium temporal and azimuthal covariant com- 
ponents of the 4- velocity: 



-£ u t . 



(22) 
= -1, 



A ddin g the normalization condition 

Eqs. (pOj) to ([22]) fully determines the fluid 4- velocity. 
Equivalent formula can be derived for all other deforma- 
tions. 



4. Light curves and power spectra of the deformed 
torus 

4.1. Ray tracing on the deformed torus 

Light curves and power spectra of the deformed torus are 
comput ed by using the gene ral relativistic ray-tracing code 
GY0T0 (Vincent et al. 20111. Null geodesies are integrated 
backward in time from a distant observer's screen to the 
optically thick torus. The zero of the surface function is 
found numerically along the integrated geo des ic and the 
fluid 4- velocity is determined by using Eqs. (20 1. The spe- 
cific intensity emitted by the undeformed torus is chosen 
to be constant throughout the surface. An image is then 
dchncd as a map of specific intensity. 

Figure [l] shows the undeformed torus image as seen from 
an inclination of 45° and 85°. The two values of inclina- 
tion emphasize the bigger impact of the beaming effect at 
high inclination, which strongly increases the dynamic of 
the image. Moreover, the apparent area of the primary and 
secondary images vary, depending on inclination. The pri- 
mary image dominates clearly at small inclination, whereas 
the secondary image's apparent area becomes important at 
high inclination. The value of f3 is chosen to ensure that 
the torus is indeed slender, i.e., that for all points inside 
the torus 



< 1, 



ro 



7T 

«2- 



(23) 



We choose /3 = 0.05, a value that satisfies these con- 
straints well to be able to still abide by the slender torus 
approximation even for the deformed (and, in particular, 
extended) torus. This value will be fixed for the remainder 
of this article. 

When the torus is deformed, two main effects will have 
an impact on the light curve: the change with time of pro- 
jected emitting area and the varying amplitude of relativis- 
tic effects as the torus surface is approaching or receding 
from the black hole. To be able to compare the respective 
impacts of these two effects on the various deformed tori, 
we require that: 



the specific intensity emitted at the surface of the op- 
tically thick torus is inversely proportional to its cross 
section. This will erase the effect of varying the flux by 
changing the torus cross-section area for torus expan- 
sion (all other deformations leave the cross-section area 
unchanged). Only the effect of changing the projected 
emitting area will thus remain; 

the 'closest approach' of the torus to the black hole, as 
defined by the minimum value of the r coordinate, is the 
same for all deformations. This ensures that the light 
emitted by the torus experiences comparable general 
relativistic effects, whatever the deformation. Precisely, 
we require that at closest approach, x —2.5 (com- 
pared with the initial minimum value of x = —1.5 for 
the undeformed torus). We note that the slender condi- 
tion in Eq. ( 23 ) is still satisfied. 



The last condition will also impose the value of the de- 
formation A parameter on all deformations. We stress that 
A is different for different deformations to allow a common 
closest approach for all deformations. 

4.2. Light curves and power spectra 
4.2.1. Definitions 

Light curves are easily computed by calculating images of 
the deformed torus at various times. The time sampling 
was chosen to be around 100 frames per Keplerian period, 
and images are 1000 x 1000 pixels. As the Keplerian period 
at radius ro = 10.8 M is tx = 230 M, the time interval 
between two frames is around St = 2 M. The total number 
of frames computed for one given light curve is N = 200, 
thus covering around two Keplerian periods. One point on 
a light curve is simply equal to the summation of all pixels 
of a given frame (this boils down to summing intensity over 
all solid angles and thus to computing a flux). 

Power spectral densities are computed in the following 
way. Let F(tk) be the flux value at time tk — kdt. The power 
spectral density at frequency fk — k/(N5t) is defined as the 
square of the modulus of the fast Fourier transform of the 
signal, thus: 



PSD{j k 



N-l 

F(t k )exp(2irijk/N) 

3=0 



(24) 



4.2.2. Observable characteristics of various deformations 



Figure |2p,) displays the light curves obtained for all defor- 
mations of the slender torus, seen with an inclination pa- 
rameter of 5°, 45°, and 85°. The panels b) and c) of the 



4 



Mazur et al.: Oscillating tori for modeling QPOs 





Fig. 1. Undeformed torus seen from an inclination of 45° (left) and 85° (right). The color bar, common to both panels, 
shows the ratio A = (t'obs/^cm) 3 , where i/ obs and v cm are the observed and emitted frequencies. This is the redshift factor 
that modulates the observed intensity. At high inclination, the beaming effect (which makes the redshift factor brighter 
on the approaching side of the torus) is more pronounced and the apparent area of the secondary image is bigger. 



same figure show the power density spectra of all defor- 
mations. The spectra are presented in two distinct panels 
for visibility as the power associated with different defor- 
mations or the same deformation at different inclinations 
varies significantly. 

The biggest modulations (or biggest power densities) 
are obtained for deformations that lead to the biggest 
change of the torus apparent area. These deformations de- 
pend on inclination, but always include expansion and one 
or more shears. 

The variation of relativistic effects linked to the chang- 
ing torus location is less important than the change of ap- 
parent area. This is clearly demonstrated by noting that the 
smallest modulation is always obtained for vertical transla- 
tion, which indeed leads to the smallest change of apparent 
area. 

Whatever the inclination, expansion leads to a high 
power density, while translations and rotation give low 
power densities. On the other hand, shears exhibit a very 
inclination-dependent power density: radial and pure shears 
give high power for low inclination, while vertical shears 
give high power for high inclination. All these observable 
characteristics are directly linked to the change in apparent 
area of the deformed torus at various inclinations. 

The time average of the light curves differs from one 
deformation to the other. It is dictated by the average value 
of the apparent area. 

The dominating frequency of the light curves depends 
strongly on the deformation. Power is balanced between 
one time and two times the Keplerian frequency i>k- 
Translations and expansion give a single-peak spectrum 
centered on v-^. Rotation also gives a single-peak spectrum, 
but centered on 2 ^k, as all possible values of apparent ar- 
eas are covered after only half a Keplerian period in this 
case only. Shears give rise to a double-peak structure in the 
power spectrum, with a different balance of power between 

and 2 vk for radial, vertical, and pure shear. 



4.3. Discussion 

Figure [2] shows that the power spectral density is an inter- 
esting probe of the deformation of a slender accretion torus 
surrounding a black hole. Different kinds of deformations 
lead to significantly different power spectra. 

It is interesting to note that the power density associ- 
ated with expansion and the dominant shears is always one 
to two orders of magnitude higher than the power density 
associated with translations and rotation. This is due to 
the highly different evolution of the torus apparent area for 
these kinds of deformations. 

The single- or double-peak nature of the power spec- 
trum is also an interesting probe of the underlying kind of 
deformation. However, the very inclination-dependent as- 
pect of the shears power spectra, as opposed to the other 
kind of deformations, presents a difficulty. 

Taking these two probes into account (maximum power, 
dominating frequency) it is possible to constrain the kind of 
deformation by studying the power spectrum. Knowing the 
maximum value of power density, it is possible to determine 
whether the underlying deformation is within {expansion, 
shear} or {translation, rotation}. Then, as a function of the 
value of the dominating frequency, it is possible to deter- 
mine the kind of deformation: expansion, shear, translation, 
or rotation. But it is not straightforward to determine what 
kind of translation or what kind of shear unless the incli- 
nation parameter is already constrained. 

However, this work is only a first step towards a more 
sophisticated treatment of the observable characteristics of 
non-equilibrium accretion slender tori. Therefore, it is not 
the aim of this article to propose a way of constraining 
the motion of the accretion torus from the power density 
spectrum properties. 

5. Conclusions and perspectives 

We compute light curves and power density spectra 
of slender optically thick accretion tori surrounding a 
Schwarzschild black hole, which are subject to simple de- 
formations (translations, rotation, expansion, and shears). 
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Fig. 2. a) Normalized light curves for torus deformations, b) Normalized power spectral densities. Set of deformations 
with large maximal normalized values, P max ^ 0.5. b) Normalized power spectral densities. Complementary set of 
deformations with small maximal normalized values, P max ^ 0.25. We note that the power spectra are represented as 
triangular functions for readability only, although they are actually made of a succession of infinitely thin peaks. 



We show that the power spectrum can be used to constrain 
the underlying kind of torus deformation. 

This is a first step towards the full study of the ob- 
servable characteristics of oscillating slender accretion tori. 
Our conviction is that combining simple deformations of 
slender tori will allow realistic oscillations to be mimicked. 
From this perspective, our present result that it is pos- 
sible to differentiate various kinds of simple deformations 
from the observed power spectra is promising. Future work 
will determine whether this still applies when combinations 



of simple deformations are considered and whether this 
could lead to a practical tool for studying power spectra 
of high-frequency QPOs to constrain the underlying physi- 
cal model. We believe that such tests can finally bring im- 
portant outputs in the field of verifying the predictions of 
general relativity in the exploration of observational data 
accumulated so far. We also believe that outputs of sophis- 
ticated ray-tracing simulations could be key components of 
the proper interpretation of data received from planned fu- 
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ture X-r ay missions such a s Large Observatory For X-ray 
Timing ( |Feroci et aL]|2012[ ). 
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